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ABSTRACT 

In this letter we explore the environment of Pluto and Charon in the far infrared 
with the main aim to identify the signs of any possible dust ring, should it exist in 
the system. Our study is based on observations performed at 70fim with the PACS 
instrument onboard the Herschel Space Observatory at 9 epochs between March 14 
and 19, 2012. The far-infrared images of the Pluto-Charon system are compared to 
those of the point spread function (PSF) reference quasar 3C454.3. The deviation 
between the observed Pluto-Charon and reference PSFs are less then la indicating 
that clear evidence for an extended dust ring around the system was not found. Our 
method is capable of detecting a hypothetical ring with a total flux of >^3.3 mjy at a 
distance of ~153 000 km (-^8.2 Pluto-Charon distances) from the system barycentre. 

We place upper limits on the total disk mass and on the column density in a reasonable 
disk configuration and analyse the hazard during the flyby of NASAs New Horizons in 
July 2015. This realistic model configuration predicts a column density of 8.7x10“^° 
g cm“^ along the path of the probe and an impactor mass of 8.7x10“® g. 

Key words, planets and satellites: rings — minor planets, asteroids: individual 
(Pluto,Charon) 


1. Introduction 


The Pluto-Charon system harbours at least four smaller satellites. Nix and Hydra 


(Weaver et al 


2006; 


Buie et al 


less, while Kerberos and Styx (Showalter et a' 


2006) haye diameters estimated to be around 100 km or 


Studies of the phase space (INaev et al 


Zsigmond fc Siili 


201C; 


Hires et al 


200i 


2011 


i 


2011 


2012) have diameters less than 40 km. 


Siili^Zsi£mpnd ^ 


Giuliatti et al 


2000 ; 


Giuliatti et al 


2010 ; 


2013h showed that dynamically 


stable regions exist around the system. Formation of a debris disk as a result of collision of 


* Herschel is an ESA space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from NASA. 
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small bodies was also suggested (lDurda^_^_Stera 


2000 ; 


Steffletal 


2006). The existence 


of such a ring was investigated bv IStefll fc: Sternl (|2007l) but observational evidence for the 
ring was not found. They placed a Scr upper limit on the azimuthally averaged normal I/F 
of ring partic l es of 5.1x10^ at a distance of 42 000 km from the Pluto-Charon barycentre. 


Boissel et al 


(|2014ll derived an upper limit of 15 000 for the number of bodies smaller than 


0.3 km at distances smaller than 70 000 km from Plu to’s system barycentr e. The impactor 


velocity of colliding bodies ranges from 10-100ms ^ (jSteffl fc Sternl . 120071) . Given that the 


escape velocities of Nix and Hydra are between 30 and 90ms“^, any material ejected from 
these bodies can escape from the satellites but cannot escape from the system (Pluto and 
Charon escape velocities are ^ 1200 and ^500ms~^, resp ectively). Recently a ring system 


around the centaurs Chariklo (iBraga-R ibas et al.L l2014^ wa s discovered and several indi¬ 
cators show that probably (2060) Chiron ( Ortiz et al.[ 2015 1 harbours a ring system, too. 


These suggest that dust rings around small bodies may be a common feature. 

The existence of a dust ring aro und the system would also mean a hazard for NASA’s 
’New Horizons’ probe (ISterrJ . 120081) . The goal of the mission is to provide the first in situ 
exploration of the system, with a close fly-by in July 2015. 


Visual wavelength range reflected light observations (IShowalter et al.l . l201lLl2012l) taken 
with the Hubble Space Telescope (HST) are, however, not ideal for the detection of diffuse 
material in the Pluto-Charon system due to the very high brightness of Pluto and Charon 
and a possibly low surface brightness of the hypothetical ring. In the thermal infrared 
the contrast between the central object (star) and the surrounding material (debris disk) 
is higher than in the visual range, due to the very different spectral energy distributions 
(SEDs) of the two components. In the case of Pluto, the thermal infrared may still be 
beneficial in the case of dark gr ains in the hypothe tical ring, like e.g. the ones found in the 
Uranus system (albedo of ~6%, 


Karkoschka. 


200111 . 


In this letter we investigate thermal emission (70/rm) images of the Pluto- 
Charon system , taken with the Photodetector Array Ca mera and Spect r omet er (PACS, 


Poglitsch et al 


20101) of the Herschel Space Observatory (IPilbratt et al.L 120101) . in order 


to identify any extended emission, originating from sources other than Pluto and Charon 
themselves. 

Diffuse dust structures like disks and rings modify the PSF of their parent sources, 
even if they cannot be fully resolved. Such partially resolved features - if bright enough 
- cause a broadening of the PSF of the source. This kind of broadening has been used to 
detect part i ally r esolved debris disks with the PACS photometer in many cases (see e.g. 


Moor et al 


2015^ . In the case of PACS, the 70p,m band has the highest spatial resolution 


(FWHM of 577), and also the SED of large grains peaks around this wavelength at the 
heliocentric distance of Pluto. In the case of smaller grains the temperature is higher than 
Pluto’s, causing a shortward shift of the peak of the SED, again favouring the shortest 
wavelength PACS band to detect it. The ratio of the fluxes from a disk/ring and Pluto is 
expected to be the highest in this band as well. Therefore we only use the TOfim band data 
in the following analysis. 
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2. Observations and data reduction 


The PACS camera on board the Herschel observed the Pluto-Charon system multiple 
times. We used the lig ht-curve measurement of Lellouch et al. (proposal ID: OT2_elellouc_2, 


Lellouch et al 


(2013) & Lellouch et al. 2015, in prep.) taken at 70/160/rm between March 


14 and 19 in 2012. All measurements were made in scan-map mode with a scan speed 
of 20" and lasted for 286 s. Pluto-Charon was observed at a heliocentric distance of 
r = 32.2 au, and at a distance of A ~32.4 au from the telescope. Basic characteristics of the 
Pluto observations are listed in Table lA.ll 

In addition to the target measurements, we used the engineering observations of the 
fiducial star a Boo for calibration purposes, as well as 42 Herschel/PACS mea surements 


Wehrle et al. 


2012) as point 


of the quasar 3C454.3 (proposal ID: TOO_awehrle_2, see 
spread function reference (discussed in details in Sect. E]). The quasar observations were 
taken between November 19, 2010 and January 10, 2011. Each observation was 220 s long, 
in the same observational mode and with the same scan-speed as in the case of the Pluto 
measurements. 

All observation s were processed with the Herschel Interactive Processing Environment 
(HIPE, lOttl . I 2 OIOII . version 13 (developer release). Starting from the Level -1 stage the 


obser vations were first corrected for the evaporator temperature effect (see 


Moor et al 


201411 . with correction factors of <0.3% to 3.2 %. Then we f ollowe d the standard pipeline 
steps to eliminate the 1 /f noise as described in 


Balog et al. 


(1201411 using iterative masked 


high-pass filtering of the observational timelines with the photProject() task. For large 
scale extended emission studies the mapmaker algorithm JScanam is often more suitable, 
but in our case the structures are kept in the inner ^30"of the source, so there is no 
need to us e JScanam to deal w ith a broadened PSF. For more details see the mapmaker 
report of (IPaladini et al.l . 12013^ . We mosaic the images with 1.1" pixel sizes. As Pluto 
was moving relatively fast w.r.t. the sky reference frame (^2.1"h“^), we constructed the 
final images in Pluto’s co-moving reference frame for each individual OBSID. We also 
used temperature correction in order to eliminate effects caused by the telescope mirror 
temperature variations. We eliminated the slight pointing jitter by using the experimental 
tool that corrects the pointing based on the information provided by the gyroscopes of 
the satellite. To alleviate the ef f ect of correlated noise we used the drizzle parameter of 

for more details on the PACS photometer noise 


pixfrac = 0.2 (see iPooesso et al 
characterisation). 


2011 


3. Thermal emission model of the Pluto-Charon system 

The observed PSF of any object is a superposition of monochromatic PSFs over the wave¬ 
length range of the bandpass filter. The shape of the object’s SED determines the observed 
PSF FWHM. Objects with high temperature have a decreasing SED in the PACS 70 fim 
band and monochromatic PSFs of the shorter wavelengths are dominant. A constant or 
rising SED at these wavelengths obviously leads to the dominance of longer wavelength 
PSFs, resulting in a broadened FWHM value with respect to e.g. that of a normal star. 
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Fig. 1. Normalised intensity of 3C454.3 subtracted from the normalised intensity of the 
object PSFs, as a function of distance from the PSF centroid. The reference profile is the 
PSF of quasar 3C454.3 (black solid [0,0] line). The calibrator star profile (blue dashed line) 
is the narrowest. The Vesta PSF (dashed magenta) is wider than the calibrator star, but 
still narrower than the reference PSF. The filled grey area represents the uncertainty of 
the 9 Pluto observations. The mean median (MM) Pluto profile is shown as the green solid 
line. The most realistic dust ring model detected at 1+3(T signihcance level is plotted as 
the red dotted line. 



In the thermal infrared the SED of Pluto is, indeed, very different from that of a 
stellar photosphere, due to the very different temperatures. Calibration stars of the PACS 
photometers have temp eratures of 3000 K or h igher, while the surface temperature of Pluto 


and Charon are ^50 K ( Lellouch et al 


201111 . Even the bright mino r planet Vesta, whose 


PSE was widely used in different Herschel applications (|Luta . l2012ll . has a typical surface 
temperature of ~200K, notably higher than that of Pluto. The spe ctral index a o f Plut o 


is ~0.82, assuming a mean Tb of 46.5K and 43K at 70 and 160 (|Lellouch et al.l . I2013ll . 
where a{iy) = For comparison, the spectral indices of a stellar photosphere and 

that of Vesta are ~2 at this wavelength. 

To exclude the possibility of misinterpreting the SED-broadened Pluto PSF as a signa¬ 
ture of a ring, we have collected the observations of the quasar 3C454. w hich is known 
to have a spectral index of 0.88 ± 0.07 at 70 - 160 /rm (|Wehrle et al.L 1201211 . very similar 
to that of Pluto. The reference PSF was created by combining the 42 observations of this 
target. The cosmological distance of 3C454.3 ensures that the region where its thermal 
emission originates from is small enough that it cannot be resolved by PACS at 70 /xm, 
hence can be used as a low temperature point source reference PSF. Using this combined 
PSF a high-resolution PSF template with 0'.'22 pixel size was created. This template is used 
in the following steps to construct the synthetic PSFs (see below). The combined quasar 
PSF has been compared with the combined PSF of the fiducial star and we found that the 
fiducial star PSF is obviously narrower (see Fig. |T] and Fig. 1C. II) . 


^ http://herschel.esac.esa.int/hpt/publicationdetailsview.do? 
bibcode=2012ApJ...758...72W 
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For each epoch a synthetic Pluto PSF was constructed using the combined quasar PSF 
as template. The synthetic Pluto PSF has two components, one corresponding to Pluto 
itself, the other describing the contribution of Charon (see Eq. [T|) . The relative positions 
of these two bodies were taken from the Horizons database, as listed in Table lA.ll and 
as presented in Fig. lA.ll In the high spatial resolution model images the two sources are 
placed in these relative positions. The relative contribution of Pluto and Charon to the 
total flux are characterised by the X ratio: 


F{x, y) = XFq{x, 2 /) + (1 - X)Fg{x -xo,y- yo) 


where F^ is the normalised synthetic quasar PSF and xq and yo are the offsets of Charon’s 
position w.r.t. Pluto. The original, high resolution image is convolved with the synthetic 
quasar PSF to obtain an ’’observable” model PSF. Due to the rotational variability of 
the surface of Pluto and Charon, the relative strength of the thermal fluxes cannot be 
determined very accurately, but their ratio is clos e to 2:1 or 3:1 depending on the model 
of their thermal properties (|Lellouch et al.l . 1201111 . In our investigation we used both ra¬ 
tios, i.e X=0.67 and X=0.75, respectively. The radial intensity profiles of the synthetic 
Pluto-bCharon PSFs agree with each other very well, indicating that the position and 
relative contribution of Charon to the total flux of the system has a minor effect on the 
intensity profile. This is mainly driven by the fact that Charon is very close to Pluto (0"6- 
0"8) compared to the size of the PSF in this band (5"7). In this respect, the difference 
observed between the radial intensity curves of the individual Pluto observations should 
mostly be attributed to measurement uncertainties. Therefore we constructed a mean mea¬ 
sured (MM) radial intensity curve of the 9 Pluto measurements, using equal weights in the 
case of all measurements (note that the integration times and other conditions were the 
same for all measurements). 

The agreement of the model and observed curves are characterised by the reduced 
X^-value: 


Xr = 


1 


N — Ndof — 1 


rrE 


{yj-yT? 

s? 


( 1 ) 


where N is the number of data points, is the degrees of freedom, yf-s are the 

normalised intensity of the j-th model at i distance and j/™ are the median observed 
normalised intensities. The Si values are the standard deviations of the observed normalised 
intensities at a specific radial distance. 


4. Intensity profile results 

When the MM curve of Pluto is compared with the mean synthetic (MS) curve of the 9 
epochs, we obtain a difference with a low level of signihcance (0.8cr) between radial distances 
0"-22". This shows that the observed Pluto PSFs are very similar to the synthetic ones 
and that we are not able to detect any extended emission feature around Pluto, i.e. there 
is no obvious evidence for the existence of a notable dust ring in the Pluto-Charon system. 


5 










G. Marton et al.: Search for signatures of dust in the Plnto-Charon system 



Fig. 2. Fo (ring excess) values as a function of distance from the Plnto-Charon barycentre. 
Red and green colours correspond to ¥piuto and Fcharon ratio of 2:1, and 3:1, for Scr 
significance. Blue and magenta colours refer to the same ratios, but for Icr. Solid, dashed 
and dashed-dotted lines represent the width of the modelled ring Ar = (y.'5, I'.'O and 1"5, 
respectively. Vertical solid black lines from left to right represent the orbital semi-major 
axis of Pluto moons Charon, Styx, Nix, Kerberos and Hydra. 


It is, however, an important question, what the total brightness of such a hypothet¬ 
ical disk/ring could be that would still be detectable with our method. To test this, we 
constructed a set of model images, in each case placing a ring in addition to the contri¬ 
butions of Pluto and Charon, as described above. The moons of the Pluto system seem 
to revolve in a well-defined plane, hence any stable ring/disk is expected in this plane 
with the highest probability (INasfv et al.l . l200ffl . We model the ring as a surface brightness 
enhancement around a specific distance (rp) from the barycentre of the Plnto-Charon sys¬ 
tem in the midplane of the moon system (tilted w.r.t. the line of sight by an orientation 
angle of 61°), assuming a Gaussian intensity profile in the surface brightness distribution. 
The parameters characterising the surface brightness of the ring are (i) rp, the distance 
of the ring from the barycentre, (ii) the width of the Gaussian Ar and (hi) F,., the total 
flux density of the ring. We constructed model images using these parameters in the range 
F^ = 0-60mJy, rp = 2.55xl0M.53xl05 km (T.T-6"6) and Ar = l.lGxlO"^, 1.74x10^ and 
2.32x10^ km (0'.'5, 0'.'75 and I'.'O). For each model we calculate the Xr of the difference. The 
Fp (ring excess) dependence of the Xr values corresponding to the Icr and 3cr significance 
levels are presented in Fig|2l where a is the standard deviation of the x^ distribution. 

Fig. m shows that our method is the most sensitive at higher radial distances, and the 
detection probability also depends on the X flux ratio. In the case of rings close to Pluto the 
detection limit is significantly higher. E.g. at the distance of Charon, we would only be able 
to detect a ring with ^100 mJy excess, therefore we are not able to constrain the presence 
/ lack of a ring within or close to the orbit of Charon. The lowest detection limit belongs 
to X = 0.67 with ring distance of 6'.'6, in this case F(lcr) = 1.8 mJy and F(3cr) = 3.3mJy. 

A possible dust ring around the Plnto-Charon system is expected to be in the re¬ 
gion where the satellites exist. We assume that the most realistic model is where X=0.67, 
rp = 2'.'2 (5.1x10^ km) and Ar = I'.'O (2.32x10^ km). In this case the ring is centered at 
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a distance between Kerberos and Nix and the ring Ar also models an extended feature 
of which 58% exists inside the orbit of the Pluto satellites. The results show that for this 
ring configuration, we obtain Fq = 6.0 and 10.8 mJy for the Itr and Scr significance levels, 
respectively. This implies that the detection probability for this geometry is ^68% for a 
ring with 6.0 mJy total flux density and 99.7% with 10.8 mJy. 


5. Discussion 


Our method presented above is able to detect infrared excess from a ring according to Fig. 
2. Any calculation that wants to derive dust surface density or impactor mass based on 
this excess has to rely on assumptions on the type, spatial and grain size distribution of 
the dust. 

Assuming that the material of our hypothetical ring is dominated by grains which are 
in thermal equilibrium with the solar radiation, we can give a rough estimate of the mass 
of this ring. Considering a black body temperature of T^ = 278.3 x (r/l[au])“^/^ = 48.9K 
I K^nnedY_fc Wyalj. l20 1 111 and an absorption coefficient of Kabs = 5.96 m^ kg“^ at 70/rm 
(|Li fc Draind . 120011) we can estimate the dust mass according to Eq. [5] below: 


^dust — 


A^E.(A) 

KabsB^{X,T) 


= 2.28 X lO’^ X 


-Ftp (A) 
mJy 


[kg], 


( 2 ) 


where A is the distance to the observer and By(Td) is the Planck function. Applying Eq. 
[2] to the most realistic case put forward in Section |4l the Scr flux limit is 10.8 mJy and we 
obtain a total disk mass of 2.5x10^° kg, equivalent to the mass of a '^SdOm-sized small 
body assuming an average density of p= 1.2gcm“^. 

The flyby of New Horizons is planned to take place inside the Charon orbit. Taking the 
maximum separation of 0"85 between Pluto and Charon we calculated the fraction of mass 
and the corresponding surface density for the area between the two objects. The fraction of 
the total ring dust mass located within 0'.'85 is 4.3x10“^ (1.1x10^kg). This corresponds to 
a surface density of 8.7x10“^° gcm“^ and an impactor mass of ^8.7xl0~^ g. This is very 
simila r, only ^10% lower, than the impactor mass of 10“^ g, obtained by 


Steffi fc Stern 


( 2007^ based on Hubble Space Telescope observations. We note that a thinner ring with 
the same brightness would result in a significantly lower impactor mass. 


6. Summary 

In this paper we searched for far-infrared signatures of a possible dust ring around the 
Pluto-Charon system. We found no clear observational evidence for a dust ring but we 
were able to put model-dependent upper limits on the possible amount of dust emission 
that could remain undetectable by our method. For a realistic ring configuration we put 
a 3fT upper flux density limit of 10.8 mJy at 70 ^m. This corresponds to a surface density 
of 8.7x10“^° gcm“^ in the region between Pluto and Charon. Our estimated impactor 
mass for the New Ho rizons flyby is very similar (~ 10% lower) to the previous results of 


Steffi &: SternI (1200711 and suggests that the satellite will not be endangered by particle 
collisions. 
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Fig. A.l. Position of Charon (red rhombs) at the time of the observations relative to Pluto 
(black plus sign). 


OBSID 

Start Time 

PApc 

[deg] 

Sep. 

n 

1342241381 

2012-03-14 03:02:01 

59.21 

0.63 

1342241382 

2012-03-14 03:07:50 

59.21 

0.63 

1342241418 

2012-03-14 19:54:46 

104.15 

0.74 

1342241419 

2012-03-14 20:00:35 

107.63 

0.71 

1342241471 

2012-03-15 12:58:39 

140.0 

0.85 

1342241472 

2012-03-15 13:04:28 

140.0 

0.85 

1342241509 

2012-03-16 06:33:59 

172.32 

0.76 

1342241510 

2012-03-16 06:39:48 

172.66 

0.80 

1342241620 

2012-03-17 00:10:12 

-145.30 

0.66 

1342241621 

2012-03-17 00:16:01 

-143.43 

0.63 

1342241655 

2012-03-17 17:31:54 

-90.0 

0.68 

1342241656 

2012-03-17 17:37:43 

-93.19 

0.65 

1342241699 

2012-03-18 11:04:14 

-52.03 

0.82 

1342241700 

2012-03-18 11:10:03 

-52.58 

0.77 

1342241865 

2012-03-19 04:42:46 

-18.95 

0.84 

1342241866 

2012-03-19 04:48:35 

-19.78 

0.80 

1342241928 

2012-03-19 20:44:31 

13.95 

0.70 

1342241929 

2012-03-19 20:50:20 

11.24 

0.70 


Table A.l. Basic properties of the Pluto-Charon Herschel observations. All observations 
lasted for 286s (single repetition). The columns are: (1) Observation identifier; (2) Date 
and start time of the observation; (3) Position angle of Charon w.r.t. Pluto; (4) Separation 
of Pluto and Charon; 
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la 


Flux ratio 

2;1 

3:1 

Ar 

O'.'5 

O'.'75 

I'.'O 

0'.'5 

0'.'75 

I'.'O 

Distance 

Flux [m 

Jy] 

Flux [m 

Jy] 

I'.T 

19.9 

14.7 

11.2 

24.0 

17.8 

13.5 

2'.'2 

7.5 

6.8 

6.0 

9.3 

8.4 

7.4 

3'.'3 

3.9 

3.7 

3.6 

4.7 

4.6 

4.4 

4'.'4 

2.5 

2.5 

2.5 

3.1 

3.1 

3.0 

5'.'5 

2.0 

2.0 

2.0 

2.4 

2.4 

2.4 

6'.'6 

1.8 

1.8 

1.9 

2.2 

2.2 

2.3 

3a 

Flux ratio 

2;1 

3:1 

Ar 

O'.'5 

O'.'75 

I'.'O 

0'.'5 

0'.'75 

I'.'O 

Distance 

Flux [m 

Jy] 

Flux [m 

Jy] 

I'.T 

35.1 

25.8 

19.6 

39.8 

29.1 

22.0 

2'.'2 

13.1 

12.0 

10.8 

14.7 

13.4 

12.1 

3'.'3 

7.1 

6.9 

6.6 

8.1 

7.9 

7.5 

4'.'4 

4.6 

4.6 

4.6 

5.3 

5.2 

5.2 

5'.'5 

3.6 

3.6 

3.7 

4.1 

4.1 

4.2 

6'.'6 

3.3 

3.4 

3.5 

3.7 

3.8 

3.9 


Table B.l. Total flux values of model rings of different configurations. We placed the 
model dust ring with a Gaussian intensity profile at different distances from with varying 
FWHM and Pluto-Charon flux ratios. 



Fig. C.l. Normalised intensity profiles as a function of distance from the PSF centroid. The 
quasar 3C454.3is shown with black solid line. The calibrator star profile is presented with 
blue dashed line. The dashed-dotted magenta line indicates the Vesta PSF. The solid green 
line shows the mean of the 9 Pluto observations (MM) while the filled grey area represents 
the measurement uncertainty. The faintest dust ring model detected at Sa significance level 
is plotted with dotted red line. 
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